Electron transfer in the nonadiabatic regime: Crossover from 
quantum-mechanical to classical behaviour 

Gunther Lang, ^ Elisabetta Paladino/'^ and Ulrich Weiss ^ 

^ Institut fiir Theoretische Physik, Universitdt Stuttgart, 70550 Stuttgart, Germany 

^ Dip. di Metodologie Fisiche e Chimiche per VIngegneria, Universitd di Catania, 

Viale A. Doria 6, 95125 Catania, Italy 

email: gunther@theo2.physik.uni-stuttgart.de, ehsa@nieso.if.ing.unict.it, 

weiss@theo2.physik.uni-stuttgart.de 

(Date: February 1, 2008) 

We study nonadiabatic electron transfer within the biased spin-boson 
model. We calculate the incoherent transfer rate in analytic form at all tem- 
peratures for a power law form of the spectral density of the solvent coupling. 
In the Ohmic case, we present the exact low temperature corrections to the 
zero temperature rate for arbitrarily large bias energies between the two redox 
sites. Both for Ohmic and non-Ohmic coupling, we give the rate in the entire 
regime extending from zero temperature, where the rate depends significantly 
on the detailed spectral behaviour, via the crossover region, up to the clas- 
sical regime. For low temperatures, the rate shows characteristic quantum 
features, in particular the shift of the rate maximum to a bias value below 
the reorganization energy, and the asymmetry of the rate around the maxi- 
mum. We study in detail the gradual extinction of the quantum features as 
temperature is increased. 



I. INTRODUCTION 

Electron transfer plays an important role in many chemical and biological reaction pro- 
cesses [|l]-§[. A prime example is the ultrafast primary electron transfer step in photosyn- 
thetic reaction centers, which is responsible for the high efficiency of the photosynthetic 
mechanism. The theory of electron transfer (ET) in chemical reactions goes back to Marcus 
and Levich . They were among the first who recognized the importance of the solvent 
environment. Since the donor and acceptor states are strongly solvated, the transfer of the 
electron entails a reorganization of the environment. This results in a free energy barrier 
separating reactants and products. The transition of the electron from the reactant to the 
product state can only take place if favorable bath fiuctuations bring reactant and product 
states into resonance. 

Electron transfer is conveniently discussed by considering the diabatic potential energy 
surfaces (PES) of each electronic state. In the classical limit, the transfer rate is determined 
by two factors. First, a Boltzmann factor with the activation energy Ea required for bath 
fiuctuations to the crossing region of the two diabatic PES. Second, an attempt frequency 
pref actor describing the probability for tunnelling once the levels are in resonance. The latter 
factor is mainly determined by the overlap between the electronic wave functions localized 



on different redox sites. For large intersite coupling, the reaction is adiabatic, whereas for 
weak intersite coupling, the reaction is nonadiabatic. In this paper, we restrict ourselves 
to nonadiabatic electron transfer, where we have fast thermal fluctuations to the vicinity of 
the crossing region followed by a fast transversal of this region. Then the intersite coupling 
A can be treated in perturbation theory, which leads in lowest order to the Golden Rule 
formula for the rate, k^ oc A"^ exp{—Ea/ k^T) . 

As the temperature is lowered to the quantum regime, nuclear tunnelling through the free 
energy barrier, which crucially depends on the detailed form of the spectral density, becomes 
effective and competes with thermal activation. We concentrate on spectral densities with 
a power law form at low frequencies, which has proven to be a useful characterization in 
many chemical and condensed matter applications. The appearance of quantum features in 
the ET rate with decreasing temperature has been repeatedly studied in the last decades 
(see Refs. |]8|-[T2| and references therein). However, a unified analytical treatment of the 



nonadiabatic case in the entire temperature regime is still missing. 

The purpose of the present work is twofold. First, we give a detailed study of the quantum 
signatures for different spectral densities. The most characteristic quantum features are the 
shift of the rate maximum to a bias value below the classical reorganization energy, and 
the asymmetry of the rate around the maximum. Second, we discuss quantitatively the 
gradual extinction of the quantum effects as the temperature is tuned from the quantum to 
the classical regime. In Section II, we introduce the model and sketch the properties of the 
classical Marcus rate. In Section III, we review the important Ohmic case, for which the deep 
quantum regime can be described in analytic form. Special attention is paid to the situation 
of a large bias energy between the donor and acceptor state for nonzero temperature. To 
our knowledge, this case has not been studied systematically until now, except for some 
numerical results in Ref. [|T3|. In the remainder, we consider the rate in steepest descent 
ranging from the sub-Ohmic to the super-Ohmic case. General features are discussed in 
Section IV. In Section V, we study the rate at T = and discuss the differences between 
the Ohmic and the non-Ohmic case. The two remaining sections are devoted to an analytic 
study of the crossover from quantum to classical behaviour. In Section VI, we concentrate 
on the case of small bias energies between the two redox sites. Finally, in Section VII, we 
consider the interesting strong bias regime, in which quantum features are most pronounced. 
We present a uniform description throughout the whole temperature range. 

II. THE MODEL 

As implied by the success of classical Marcus theory, the electron can be described in 
terms of two discrete states at all temperatures of interest. Modelling the solvent by a linearly 
responding heat bath, the proper Hamiltonian is the spin boson Hamiltonian [14-T6[ 



H = -|(Aa. - ea,) + J2[^ + l^au^lxl] - ^i£a, . (1) 

a 

Here, az and a^ are Pauli spin matrices. Throughout, we put e > 0. The donor and 
acceptor states are the eigenstates of a^ with eigenvalues +1 and —1, respectively, and they 
differ by a bias energy he. Transitions between these localized states are induced by the 



intersite or electronic coupling ^A. The spin operator cr^ is coupled linearly to the collective 
bath mode ii£ = {qo/2)'^^CaXa, where £ can be considered as a fluctuating dynamical 
polarization which vanishes on average, and 2/i is the difference in the dipole moments of 
the two electronic states. The effects of the solvent are fully described by the spectral 
density G{uj) = [ql/2h) 'Yj^{c^^/maOJa)5{ijJ — Ua). Outer-sphere ET reactions are controlled 
by fluctuations in the solvent polarization which cover a wide frequency range. To model 
the solvent environment, we choose a smooth spectral density with a power law form at low 
frequencies and an exponential cutoff at high frequencies. 
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Here Uc is a cutoff frequency for the bath modes and 6s is a dimensionless coupling parameter. 
For convenience, we discriminate the reference frequency cUph from the cutoff frequency Uc- 
Ohmic or frequency- independent damping is described by s = 1, whereas the regimes s > 1 
and s < 1 are referred to as super-Ohmic and sub-Ohmic damping. 

For sufficiently high temperatures and/or strong damping, the system exhibits over- 
damped exponential relaxation towards the equilibrium state. Hence the dynamics in this 
regime is completely described by transition rates which obey the principle of detailed bal- 
ance. In the nonadiabatic regime A -C tUc, the incoherent tunnelling rate from the donor to 
the acceptor state (forward rate) is given by the Golden Rule expression [|T5| , 
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(3) 



In Eq. @, the system-environment coupling is considered in all orders, whereas the elec- 
tronic coupling is treated in second order. The quantity Q{z) is the twice integrated bath 
correlation function, d?Q{z)/dz^ = {2ijL/h)'^{£{z)£{0)), and is given by UlSlJlT 
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(4) 



where /3 = Xjk^T . The function Qi^z^ is an analytic function of the complex time z in the 
strip > Im^; > —7i(5 and obeys the symmetry relation Q{—z — ih(3) = Q{z). Therefore, we 
can arbitrarily distort the integration contour in Eq. (^) inside this strip. We may choose, 
e.g., the contour C^ with constant negative imaginary part < r < hp, as shown in Fig. |[ 
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FIG. 1. The complex time plain z with three possible contours C^ 



If we had derived the rate formula (|^) with the alternative Im F- method instead of the 
Golden Rule approach, we would have directly obtained the contour Cf where f is the width 
of the bounce (kink anti-kink pair) in imaginary time. The equivalence of these seemingly 
different approaches follows from the analytic properties of Q{z). Putting r = h(3/2 and 



z = t — ifi(3/2, the rate (H) is equivalently expressed as |T^ 
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The function X{t) is real and symmetric in t 

X{t) =Q{z = t- ihl3/2) = 



°°,.m 



u^ 



coth.{^hi3uj) — 



COs(co't) 



(6) 



In the representation (^), the detailed balance relation between the forward rate k~^{e) and 
the backward rate k~{e) = k~^{—e) is directly visible, 

k+{e) = e^^'k-{e) . (7) 

With the form (§), the integrations in Eqs. (1) and (P) can be performed exactly, yielding 

Q{z) = 25sT{s - 1) f^\ U^ " (^ + ^^c^)'-') + 2{h(3u,y-%is - 1, 1 + «:) 



X{t) = 26sT{s-l)(^] 
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Here, C,{yiq) is Riemann's generalized zeta function |T8|, T{y) is Euler's gamma function, 
and K = l/%f3uOc. For integer s, Eqs. (P) and (^) can be expressed in terms of polygamma 
functions. We shall refer to the case k -C 1 as the quantum limit, whereas the case k ^ 1 
describes the classical limit in which k-^T is the largest energy scale of the problem. Note 
that for a high cutoff cjc, the quantum regime may extend up to quite high temperatures. 
Oppositely, for a very low cutoff, the classical regime reaches down to low temperatures. 

In the classical limit k ^ 1, the integral in Eq. (^ is dominated by the short-time regime, 
in which we can substitute the classical form Qc\{z) for Eq. (^, 



Qd{z) = iAciz + {Ad/h/3)z'^ . 
Here, ^A^ is the classical reorganization energy of the solvent, where 

G{u) 



A 



cl 



duj 



UJ 



26sT{s){uJc/^phy ^c ■ 



(10) 



(11) 



The second equality holds for the spectral density (H). With Eq. (|IOD , the calculation 
of the rate (^) reduces to a Gaussian integration, giving the classical Marcus rate in the 
nonadiabatic limit 



The influence of every bath is thus characterized by only one parameter, the classical reor- 
ganization energy, and the rate is insensitive to details of the spectral density. If the bias 
equals the reorganization frequency, e = e* = Ad, the transfer is activationless, and the rate 
is maximal for fixed T. Furthermore, the rate is symmetric around the maximum. 

As the temperature is lowered to the quantum regime, nuclear tunnelling through the 
free energy barrier becomes effective. The quantum rate then shows some characteristic 
features which sensitively depend on the detailed form of the spectral density G{uj). In the 
next section, we review the Ohmic case s = 1, in which the integral in Eq. (|^) or @ can be 
calculated in analytic form in the quantum regime. For non-Ohmic spectral densities, this 
is no longer possible, and we have to employ the method of steepest-descent. 



III. THE OHMIC RATE IN THE QUANTUM REGIME 

In the Ohmic case s = 1, the spectral density (0) takes the form G(uj) = 2K u e~'^^'^'', 
where the dimensionless coupling strength Si = K is often called Kondo parameter. By 
carefully performing the limit s — >■ 1 of the expressions (|^) and (^, we obtain 

Qi.) . 2K HI + >...) + 2K ,n (^^^^—^10^01—^) , (13) 

X(t) = 2K In (-, ^l^,i;(lpl _) . (14) 



In the quantum regime h/3u)i. 3> 1, the second term in Eq, (|T^ is simplified and we get 

Q(.) = 2A- ln(l + ,.,z) + 2K hr (^HhMm) . (15) 

\ nz/np ) 

The first term represents the correct T = limit without any approximation. In the often 



employed scaling limit, the first term in Eq. (|T^) is approximated by 1K\\\.{iiOf.z). The 
scaling form for Qi^z) is only justified if the contribution of the core region \z\ < l/uc to 
the rate integral (0) is negligible. This is the case only if the bias is much smaller than the 
cutoff frequency, e ^ c^c. The equivalent scaling form of X{t) reads 



71 \npj 

With Eq. (|16|), the rate (|^) is calculated in analytic form as 

,2 /fc,Q, , \ 1-2X 



(16) 
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which is vahd for k-gT, he -C huc- The behaviour of this aheady well-known expression is 
discussed in Refs. |15,16]. We remark that the form ( p!7|) was originally obtained from the 
representation (|^) using the scaling form of Q{z) [p!9|-[2T|]. In this representation, however, 
the rate integral (|^) seems to exist only for K < ^, since the integrand is singular at the 
origin. Actually, this singularity is spurious due to the above analytic properties of Q{z), 
and therefore Eq. ([17D is valid for all K > 0. In the representation (1) with Eq. (p!6D , the 
singularity is far from the integration contour for all K. 

So far in the literature, little attention has been paid to the case where the bias is of 
the order of Uc or larger, and the bath is still quantum mechanical, h(3uJc ^ 1 |]T3[. If the 
bias is of the same order of magnitude as Uc, most of the contribution to the integral (^ 
comes from the core region \z\ < l/cjc- Then the scaling form of Q{z) is inappropriate and 
we have to use the expression (pISl). With this form, however, it is not possible to express 



the integral in (^ in analytic form for arbitrary temperatures. At zero-temperature, where 
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can be performed exactly, yielding [pl| , p2| 
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The core region 1^1 < l/cj^ gives rise to the factor e~'^/^''. For e <^ c^c, this factor is unity 
and we recover the T = limit of Eq. (|17D . For a larger bias, the factor e~'^/^'' represents 
an essential modification of the T = limit of the scaling result (|1^. For K > 1/2, the 
rate has a maximum when the bias equals the value e* = {2K — 1)uJc = Ad — ujc- In the 
second form, we have used the reorganization frequency ( ]TT| ) in the Ohmic case. Ad = 2KuJc. 
Thus the position of the rate maximum is at a lower bias value compared to the classical 
regime, where e* = Ad. The shift of the maximum is due to quantum mechanical nuclear 
tunnelling through the free energy barrier, and the relative shift c^c/Ad = 1/2K is small for 
large coupling strength K. For K < |, no rate maximum exists. However, depending on 
the magnitude of the bias, electron tunnelling might be coherent and thus a rate description 
may be inappropriate in this coupling regime. 

At finite temperature, the contribution of the core region to the integral (j^) can not be 
given in analytic form. However, the exact low temperature expansion for arbitrary bias 
can be found. We proceed by writing Eq. ( ]T5| ) as Q{z) = 2K ln(l + iojcz) + 5Q{z). In the 
regime K k^T <^ he, we may expand 5Q{z) into a power series of {z/h/Sy, 
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Substituting the series (p!9| ) into the integral (|]), we see that every term of the resulting series 
for exp[—6Q{z)] can be generated by suitable differentiations with respect to the bias. Thus, 
in principle, the full low temperature expansion can be constructed upon using Eq. (|TBp. In 
leading order, we obtain the exponential correction 
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(21) 



which varies as exp(aT^). Depending on the parameters, the sign of a can be positive or 
negative. The contribution 1/uJc in the round bracket conies from the core region. This term 
is neghgible when the bias is considerably below {2K — 1)uJc-, and in this regime Eq. ( pT]) 
coincides with the low temperature expansion of the rate (p!?!). For a stronger bias, however, 
this term is of crucial importance. The maximum of the rate (pT]) for i^ > | is at e = e* 
where 



e* = {2K - l)cjc + 
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With the onset of thermal activation at finite T, the rate maximum is shifted to a higher 
bias. 

Let us now expand the rate expression (^) around the maximum, Eq. (^21). In contrast 
to the classical formula (|T^, we now have an asymmetric distribution 
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The dimensionless width W and the asymmetry parameter A are given by 



W^ = {2K-l){l + 
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(24) 
(25) 



For K ^ 1, Eq. (|23|) is accurate over a wide bias range. Since A > for i^ > i, the wing of 
the rate above the maximum is enhanced compared to the wing below the maximum. This 
is due to the smaller width of the free energy barrier in the inverted regime e > e*. At T = 0, 
we have W^ = {2K — 1) and A = 2/[3{2K — 1)]. The width grows with increasing damping 
strength K, whereas the asymmetry decreases. For finite temperatures, thermal activation 
over the free energy barrier becomes possible. For K > 1, this leads to a widening of the 
distribution and a reduction of the asymmetry. The quantum features are reduced and the 
system behaves more "classically" with increasing temperature. In the regime K < 1, the 
behaviour is opposite: the width decreases and the asymmetry grows for finite temperatures. 
However, for K < 1 the representation ( ^31) is only accurate in a small interval around e*. 

So far we have studied the Ohmic case in the quantum regime k -C 1. It is now interesting 
to see whether the quantum characteristics change qualitatively for a non-Ohmic spectral 
density. 



IV. EVALUATION OF THE RATE IN STEEPEST DESCENT 

For a general spectral density, it is not possible to calculate the rate (0) or (||) in analytic 
form. To evaluate the integral asymptotically, we choose the contour Cf which passes through 
the saddle point Zs = —if in the direction of steepest descent (cf. the discussion in Section 
II and Fig. I). Putting 



Q{z) — iez 



(26) 



the saddle point Zs is determined by the transcendental equation 

F\z,) = , (27) 

where the prime denotes differentiation with respect to the argument. Now, expanding F{z) 
around Zg up to second order, 

F{z) = F{zs) + '^F"{zs){z - Zsf + 0[{z - zs)% (28) 

we arrive at the steepest descent expression for the rate 0,|lll, 

AV 27r ^"' 



k-'iT, e) = — [jJi^^) exp[-F(z.)] . (29) 

The quantities F{zs) and F"{zs) are positive real. Here, (3~^F{zs) is the "quantum activation 
free energy" of the donor state, and F"{zs) introduces quantum effects in the pre-exponential 
factor. As we shall find, the factor F"{zs) is responsible for the shift of the rate maximum to 
a bias value e* below Ad. The leading correction to the formula (29) comes from the fourth 



order contribution in Eq. (^). This contribution is small when the condition 

Q^^(z.)/8[Q"(^,)]2 « 1 (30) 

is fulfilled. Physically, this is the regime of multi-phonon emission processes. 

In the absence of a bias, the transcendental equation (^) reduces to Q'{zs) = 0. With 
the general form (^, the stationary point is then given by Zs = —ifij3/2. As the bias is 
increased, the stationary point moves from Zg = —ih(3/2 towards 2;^ = 0, and in general 
one has to solve the transcendental equation (P7|). The stationary point finally reaches the 
origin if the bias equals the classical reorganization frequency, e = A^ = —iQ'ifi). 

In the classical limit h(3ujc <^ 1, the stationary point is always near to the origin. Then 
the expansion (p!OD around z = is justified for any bias. 

V. ZERO TEMPERATURE RATE 

In this section, we discuss the modifications of the quantum rate features when the 
spectral density is changed from Ohmic to sub-Ohmic and super-Ohmic. 
For T = 0, the pair interaction (|^) reduces to the form 

Q{z) = 26sT{s-l){ujJuo^^y-\l-{l + iuj,zf''']. (31) 

Substituting Eq. (0) into Eq. (|26|) , the saddle point is found in analytic form, Zg = 
— «[ (Aci/e)^'* — l]/ujc. Then the rate is obtained from Eq. (^) as 
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The condition (|30D for the validity of the form ( p2D reads 

[25j:{s) ]i/^(e/cuph)(^-^)/^ = (e/Aei)(^-i)/^AdM » {s + 2){s + l)/8s . (33) 
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This roughly corresponds to the regime e < cjph in the sub-Ohmic case, e > tUph in the 
super-Ohmic case, and to 5^ = if S> 1 in the Ohmic case. It can easily be seen that the 
limit s — >^ 1 in Eq. ( |32D reproduces the Ohmic rate (p!8| ) for K ':^ 1. The rate ( |32D has 
a maximum when the bias takes the value e*, which is determined by the transcendental 
equation 
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(34) 



The expressions (0) - (|34D hold for general s. 

In the sub-Ohmic case s < 1, the formula ( P^ is valid down to e = 0. We have 
k^(e — > 0) oc e"*^**"''"'^-*/^* exp [—as(Aci/e)^-^~^^^^]. As a result of the high density of low-energy 
excitations, the barrier crossing is completely quenched for a symmetric system at zero 
temperature. In the limit Uc —>■ oo, the position of the rate maximum is selfconsistently 
determined from Eq. (0) as 



e* = [2s/(l + s) 



ls/il~s) 



[25sTis) 



11/(1-^) 



LU- 



ph 



(35) 



In the descending flank of k~^{e) above e*, the multi-phonon emission processes gradually 
die out as e is further increased. Then the steepest descent rate (|3^) ceases to be valid. For 
e ^ e*, the influence of the environment is weak and emission of energy is determined by 



the one-phonon process |16 



The super-Ohmic case shows opposite behaviour. Here the low bias rate is determined 
by one- or few-phonon processes. With increasing e, multi-phonon processes come into 
operation and form the ascending flank of A;"'"(e). For e > tUph, the asymptotic multi-phonon 
expression (|32D is valid. For Ad ^ Uc, the position of the rate maximum is found from 
Eq. (HD as 
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which correctly reproduces the Ohmic result |^. For fixed cUc/Ad, the value of e*/Aci slopes 
down with increasing s. Expanding the rate around its maximum e* as in Eq. (^3]) , we 
obtain for the width W and the asymmetry parameter A 
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In the limit s — > 1, the width (^) reduces to the Ohmic form (p^, and the asymmetry 
parameter is consistent with Eq. (|25|) for i^ ^ 1. Interestingly, the width grows and the 
asymmetry decreases with increasing s for fixed cUc/Ad. 



VI. QUANTUM TO CLASSICAL CROSSOVER FOR SMALL BIAS 



The ET rate for a symmetric system (e = 0) has been the subject of several works using 
different analytical approximations and numerical methods [p|-[TTl,p^,p^ . Here we study 



the rate in the unbiased case for a general spectral density employing the steepest descent 
method. We are particularly interested in the crossover from the quantum to the classical 
regime. Moreover, we show how the procedure can be extended to the case of a small bias. 
Since the stationary point is at Zg = —ih(3/2 for zero bias, it is natural to switch to the 
representation (^, where the integrand is stationary at tg = 0. Expanding the integrand of 
the representation (||) about t = 0, we find 

X{t) = (/i/3/4)Ai + (A2//1/3) f + 0{t^) , (39) 

where the frequency scales Ai and A2 are given by 



n(3 Jq uj^ 



Ai = — / du — ^ tanh{\h(3u) , (40) 



A, = ^ TduGiu) . I (41) 
/ Jo smh[^n/JU!) 

For the spectral density (H), we obtain the analytical expressions 

Ai = SBsik^T/h) + 6,(/c) (kBT/h) {T/T^^y-' , (42) 

A2 = as{K) (kBT/n) {T/T^^y-^ , (43) 

6,(ft:) = mj:{s - 1)|C(S - 1, 1 + «:) - C(s - 1, I + ft:)} , (44) 

a,(fi:) = 25,r(s + l)C(s + l, i + K), (45) 



where Tph = hupi^/k^ and B^ = 6sT{s — l)(u;c/^ph)^"^- There follows from the forms (ITll), 
(^), and ( PD that A^ > Ai > A2. In the limit k — > 00, the energies ^Ai and hA2 approach 
the classical reorganization energy, 

lim Ai = Aci , lim A2 = Ad . (46) 

K— »00 K— >00 

The limits (|^) are found either from the integral representations or by substituting the 
asymptotic expansion of the zeta function into Eqs. (|4^) and (143|) . In the quantum regime 
K <^ 1, the coefficient functions take the simpler form 

6,(0) = 166sT{s - 1)[2 - 2^-i]C(s - 1) , (47) 

a,(0) = 26sT{s + lyr^' - l]C{s + 1) . (48) 



Substituting Eq. (^) into Eq. (|^), the rate for zero bias is obtained as p|, [l0|Jl4 



k^T,e = 0) = ^^e-^P^^/\ (49) 

This formula is a quantum generalization of the classical result (|l^) for e = 0. The quantum 
features are included in the effective solvation energies hAi and /1A2. The energy hAi/A plays 
the role of a quantum activation free energy, which is below the classical value Ea = /iAci/4 
because of nuclear tunnelling. The energy ^A2 accounts for quantum effects in the attempt 

10 



frequency due to bath fluctuations. Both quantum effects lead to an enhancement of the 
rate over the classical value. The leading correction to the formula (^) arising from the 
(9(t^)-term in the expansion (^) is small when the condition 

is satisfled. Physically, this condition holds in the regime of multi-phonon processes. It is 
always fulflUed in the classical limit. In the quantum regime k <^ 1, Eq. (|50|) corresponds 
in the Ohmic limit to the damping regime K ':^ 1. In the super-Ohmic case, the condition 
(0) is met for temperatures well above Tph, whereas for a sub-Ohmic bath, it is fulflUed 
when T is fairly below Tph. 

For flnite bias, the stationary point moves away from t^ = 0. However, for a small bias 

e<(l + 2fi:)A2, (51) 

the short-time expansion (^) for X(t) can still be applied. The stationary point is located 
at ts = ihPe/ 2 A2, and the rate (|) is found as 



*^(T.^)-^,^e.p(-5£Ai,5^_M£!.. (s„ 



4 V A2 i 4 2 4A2 

This expression describes the crossover from the classical to the quantum regime for general 



spectral densities. Corrections to the formula (^) are small if the conditions ( [50| ) and ( |5TD 
are fulflUed. In the classical regime k ^ 1, the condition (|5T|) is satisfled even when the bias 
is very large. We see from Eq. (|4^) that the expression (^) reduces to the classical rate (|12D 
as K — »• (X). As the temperature is lowered to the quantum regime, the particular properties of 
the spectral density G{u) become relevant which indicates that nuclear tunnelling becomes 
effective. The quantum effects are captured by the deviations of the effective solvation 
energies /lAi and ^A2 from hA^. However, the bias range ( ^1]) in which the rate formula 
(^) is valid shrinks as temperature is lowered. In the quantum regime k <C 1, the condition 
(pl|) reduces to e <^ A2 oc T"*. Hence the bias range goes to zero as T — >• 0. In the quantum 
regime k -C 1, the formula ( ^2] ) does not properly describe the rate around the maximum. 

In the Ohmic case s ^ 1, we obtain from Eqs. ( ^21) and (|43|) in the quantum regime 
K < 1 

Ai = {8K/h(3) \n{nf3uJc/7T) , A2 = n^K/n(3 , (53) 

and the rate expression ( |52D takes the form 

*^(T..^^,/^(^V"^xp|Ml-lfMiV\ (54) 



4cJc V K \fiP^cJ 12 K \2n 

Now, the conditions ( pOD and ( [5T| ) reduce to i^ ^ 1 and he <^ tt'^K k^T. In this parameter 
regime, the steepest descent-rate ( ^41) coincides with the previous expression (|l^). 

For a larger bias than given by Eq. ([5T|), the stationary point is no longer near to t = 
and the expansion (^) cannot be employed any more. This case is studied in the next 
section. 
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VII. QUANTUM TO CLASSICAL CROSSOVER FOR STRONG BIAS 



As the bias is further increased, the stationary point Zs approaches the origin of the 
complex time plain, and for e = Ad, we have Zg = 0. In general, for large bias and/or high 
temperature, the stationary point is near z = 0. To explore this regime in which quantum 
signatures are most pronounced, we expand the integral representation (^ of Q{z) around 
z = up to fourth order in z, 



Q{z) = iAciz + 



C2(l^ 



iUJcZ 



C3 



[UJrZ] - 



C4(k) 



The dimensionless coefficient functions Cj{k,) are given by 

1 



UcZ 



+ Oiu.zf 



(55) 



C2(KI 



UJ^ 



^3 = ~ 






dujujG{uj) , 

duj uj'^G{uj) coth(^hPuj) . 



(56) 



(57) 



For the spectral density @, we readily obtain 

C2in) = s[l + 2k'~^\{s + 1, 1 + k) ] (A,iM) , 

C3 = s{s + l)(Aci/u;c) , 
Ci{K) = s{s + l){s + 2)[l + 2K'+'^Cis + 3,1 + k)] (Acl/uJc) ■ 

In earlier works, the analysis of the rate has been based on truncating the series (|55D 
at order z"^ P,p!2|. However, this crude approximation, which has been often termed "semi- 
classical" , fails to describe the actual shift of the rate maximum e* away from A^ and the 
asymmetric distribution around the maximum in the quantum regime k <C 1, as noticed 



already in Ref. |2J|. In the sequel, we study the quantum effects of the rate related to terms 
of order z^ and higher in Eq. (^) . 
In the bias regime 



the stationary point is obtained from Eq. 
Up to second order, we find 



z.. 



cl 



Then the rate is found in the form (T2 



A, 



(58) 
with Eq. (|55|) in the form of a series in A^ — e. 

(59) 



J_C3 
2uJc C2 



A 



cl 



C2UJc 



A2 



Fiz^ = — 



2tx 



2co 



Aci-e 



1/2 



UJr 



+ 



exp[-F(z^)] 
Aci- 



_C3_ 

6co 



LOr 



(60) 



F'\zg) 



C2 - 



C3 
C2 



Ah - 



UJr 



C4 



A 



cl - 

OJr. 



UJr 
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This rate expression together with Eq. ([56|) describes the full crossover from quantum- 
mechanical to classical behaviour for a general spectral density. It is valid if the bias condition 
(^) and the steepest descent condition (|30|), which approximately takes the form C4/8C2 <C 1, 
are fulfilled. For /t ^ 1, we see from Eq. (|56D that these conditions are fulfilled for any bias 
and spectral density, and the rate reduces to the classical result (p!^). For lower k, the rate 
expression (|60|) shows substantial deviations from the classical behaviour. The quantum 
activation factor exp[—F{zs)] of the rate alone would predict activationless transfer for 
e = Aci. However, the rate maximum is actually shifted by the fluctuation factor F"{zs). 
We consistently obtain 



Aci 2c2Aci ' 

The rate around the maximum is again described by the form 
asymmetry parameter A are given by 



w 



A 



C2 



1 



„3 

'^2 



C4 

^2 



_C3_ 

3c2 



(61) 

^). The width W and the 

(62) 
(63) 



We emphasize that the expressions (pT| ) - ( |63D together with Eq. (|56| ) are valid for a general 



spectral density. In the following, we will restrict ourselves to the special form (^) of the 



coefficient functions Cj. Before discussing the quantum signatures in the parameters e* 



W, and A in detail, we now consider the conditions (|58|) and (pO]) in the quantum regime 
K <^ 1. First, the steepest descent condition ( ^0|) is approximately given by Ad ^ Uc- This 
condition is practically always satisfied in the super-Ohmic case when Uc ^ cUph, whereas 
in the Ohmic case, it corresponds to the regime ii' ^ 1. It is usually not fulfilled in the 
sub-Ohmic case, except for very strong coupling 26sT{s) ^ (cUc/cjph)^"*. Secondly, the bias 
condition (^) takes roughly the form |Aci — e| ^ Ad, and thus at T = we are confined 
to a narrow bias interval around Ad ^^. Importantly, the most interesting area around the 
rate maximum falls into this interval. It is easily verified that in this bias interval the T = 
limit of the rate (pDI) reduces to the previous result 



e*/Aci 



1 




0.98 


^-^TTTT^^"^^^'- 




/ ^ _,-■'" 




/ ^ ■'' 




1 ^ •' 


0.96 


1 ' 




1 ' 




/ 


0.94 


/ 
/ 




/ 




/ 


0.92 




0.9 
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FIG. 2. e*/Ad as a function of k, for aJc/Ad = 0.05. The solid, dashed, and dotted lines corre- 
spond to the cases s = 1, 2, and 3, respectively. 
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The position of the rate maximum, Eq. (|61|), is plotted as a function of k for different 
values of s in Fig. ^. In the asymptotic low temperature regime, we find from Eq. (61) 



e* s + 1 ujc 

Aci ~ 2 Ael 



i-2C(s + i: 






s+1 



(64) 



which correctly reproduces the T = result, Eq. (p^). The leading low-temperature correc- 
tions vary as T*"*"^. For s = 1, the form (|6^ reduces to the earlier result (22) for i^ ^ 1. 



With increasing k, the maximum e* is shifted upwards until it finally reaches the classical 
value Aci. With increasing s, the approach to the classical limit is more gradual. As visible in 
Fig. ^ the deviations from the classical behaviour are more distinct for larger s throughout 
the whole temperature range. 

The width and the asymmetry parameter of the rate are displayed in Figs. § and ^. 



W^ 




FIG. 3. The squared width W^ as a function of k for uJc/Ad = 0.1. The solid, dashed, and 
dotted lines correspond to the cases s = 1, 2, and 3, respectively. The thin solid line represents 
the classical component W^i = 2(Aci/wc)«^- The difference between this straight line and the full 
squared width is the quantum component W^^. 

At very low temperatures, the expressions 

sAcl / S + 1 UJc 



w 



A 



UJr 



1 



2s Kr 



3sA 



cl 



K{s + 1) 



l + 2C(s + l)(l 



and (|63D take the form 



kWT 



S — 1 CiJc 

2s A~ 



hUr 



S+1 



huJr 



s+l 



(65) 
(66) 



Again, the T = results (|37D and (|38D are correctly reproduced p6|. With increasing s, the 
width grows, whereas the asymmetry decreases. The low temperature corrections vary as 
T'^'^^. For s — > 1, we recover the Ohmic results ( p^ and ( pSf ) ioi K ^ 1. With increasing 
temperature, the asymmetry decreases and the width grows. However, the increase of the 
width is due to the "classical component" VFj = 2{Kc\/^c)i^ arising from thermal activation, 
whereas the "quantum component" W^^ decreases. For this reason, the quantum effects are 
reduced with increasing temperature. In the limit k — > cxd, the asymmetry finally drops down 
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to zero, and the width approaches the classical form W^i = 2(Aci/ct;c)/t. With increasing s, the 
curves in Figs. |^ and ^ approach the classical limits more gradually, and the deviations from 
the classical behaviour are more pronounced. This is in correspondence with the behaviour 
of e*. Thus, with growing s, the quantum regime extends to higher temperatures. Clearly, 
this is due to the larger contribution of high frequency modes leading to an enhancement of 
nuclear tunnelling. 
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FIG. 4. The asymmetry parameter ^ as a function of k for cOc/Ad = 0.1. The solid, dashed, 
and dotted lines correspond to the cases s = 1, 2, and 3, respectively. Since A > 0, the wing of the 
rate above the maximum is enhanced compared to the wing below the maximum. In the quantum 
regime k <C 1, the deviations from the classical value A = are more distinct for smaller s. This 
is in contrast to the behaviour of e* and W'^ for k <^ 1 displayed in Figs, pi and H. 

VIII. CONCLUSIONS 

We have studied the nonadiabatic tunnelling rate in the biased spin-boson model within 
the Golden Rule formula. We have analyzed the quantum signatures of the rate in the 
position of the maximum e*, the width W, and the asymmetry parameter A for different 
spectral densities. In the Ohmic case, we have calculated the rate in analytic form at T = 0, 
and we have presented the leading low temperature corrections. Since the rate maximum 
is located in the strong bias regime, it has been essential to treat the bath correlation 
function beyond the scaling limit. For an Ohmic and non-Ohmic spectral density, we have 
investigated the entire quantum to classical crossover of the rate within steepest descent, 
and we have discussed the parameter bounds within which this approximation is justified. 
We wish to emphasize that the characteristic quantum features found are not artefacts 
of the steepest descent method. For the parameter regions displayed in Figs. @ - @, the 
deviations of the analytical expressions from numerical computations of the Golden Rule 
rate are negligibly small. The quantum signatures of the rate are gradually extincted with 
increasing temperature. We found that for fixed low temperature, k <^ 1, the quantum 
features in the position of the maximum and the width are growing with increasing s. In 
contrast, the quantum signatures in the asymmetry parameter become less pronounced if we 
switch from Ohmic to super-Ohmic. At higher temperatures, k > 1, the quantum signatures 
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in the parameters e*, W, and A are more distinct for larger s. Thus in the super- Ohmic 
case, the regime in which deviations from the classical behaviour are significant extends to 
higher temperatures than in the Ohmic or sub- Ohmic case. 

Many of our results are not limited to the particular form (j^) of the spectral density. 
They also apply to spectral densities with a peak-like structure at higher frequencies, which 
are relevant when inner-sphere reorganization or specific vibrational modes in the solvent 
are important. The quantum effects studied here are particularly pronounced for solvent 
modes with energy above the thermal energy. Finally, our study is also relevant for single 
charge transfer in ultrasmall junctions in which the electromagnetic environment plays the 



role of the solvent pTlp^ 
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